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Abstract. Steady axisymmetric outflows originating at the hot coronal magnetosphere of a Schwarzschild black 
hole and surrounding accretion disk are studied in the framework of general relativistic magnetohydrodynamics 
(GRMHD). The assumption of meridional self-similarity is adopted for the construction of semi-analytical solutions 
of the GRMHD equations describing outflows close to the polar axis. In addition, it is assumed that relativistic 
effects related to the rotation of the black hole and the plasma are negligible compared to the gravitational 
and other energetic terms. The constructed model allows us to extend previous MHD studies for coronal winds 
from young stars to spine jets from Active Galactic Nuclei surrounded by disk-driven outflows. The outflows are 
thermally driven and magnetically or thermally collimated. The collimation depends critically on an energetic 
integral measuring the efficiency of the magnetic rotator, similarly to the non relativistic case. It is also shown that 
relativistic effects affect quantitatively the depth of the gravitational well and the coronal temperature distribution 
in the launching region of the outflow. Similarly to previous analytical and numerical studies, relativistic effects 
tend to increase the efficiency of the thermal driving but reduce the effect of magnetic self-collimation. 
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1. Introduction 

The formation of relativistic jets around compact ob- 
jects and Active Galactic Nuclei (AGNs) is one of the 
most intriguing and yet not fully understood astrophysi 



cal phenomena ((Ferrari 19981 Mirabel fc Rodriguez 1998 



IMirabel 2003|l . In those jets, velocities reach a fraction of 
the speed of light with the corresponding Lorentz factor 
ranging from values 7 ^ 2 — 10 in Seyfert Galaxies and ra- 
dio loud AGNs (Piner et al. 2003, |Urry fc Padovani 19951 ) 
up to the inferred values 7 = 10 3 in GRBs; AGN jets are 
also characterized by the rather narrow opening angles of 
a few degrees (Biretta et al. 2002, Tsinganos & Bogovalov 
2005). 

MHD models for coronal or disk-jets rely on the basic 
idea that the gravitational energy of the central object is 
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transferred to the accreting plasma which via a collimation 
mechanism then produces the jet. This energy released by 
accretion increases with the mass of the central object, 
a fact which may explain the wide variety of the power- 
ful jets observed. Several analytical and numerical efforts 
have been invested to investigate the mechanisms of jet 
acceleration and collimation. The formation of collimated 
jets seems to be closely related to the presence of large 
scale magnetic fields (e.g., IGabuzda 2003(1 and the exis- 
tence of a gaseous disk and/or a hot corona around the 
central object flKonigl fc Pudritz 2000"|lLivio 2002)l . 

Regarding the energy source of jets, it is usu- 
ally assumed that at their base they are powered ei- 
ther by a spinning black hole (Blandford & Znajek 
1977. IRees et al. 19821 TBegelman et al. 1984> , or, by the 
surrounding accretion disk (Miller & Stone 2000). 
Furthermore, they are plausibly Poynting flux dominated 
l|Sikora et al. 20 05 ) with their central spine hydrodynam- 
ically dominated l|Meliani et al. 2004)) and their plasma 
composed by protons-electrons or by electron-positron 
pairs. 

Regarding the collimation of the outflow, this is 
likely to be due mainly to the hoop stress resulting 
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from the toroidal magnetic field generated by the ro- 
tation of the source (Bogovalov 1995). Recent VLBI 
observations suggest that the direction of the mag- 
netic field vectors is transverse to the jet axis. This 
is the case in BL Lac objects, e.g., in Mrk 501 
uzda 20 03). or, in quasars where the central faster 
part of the jet is characterized by toroidal magnetic 
fields (Asad a~et al. 2002}) . Magnetic self-collimation has 
been shown to be efficient in the non relativistic 
context ( jHeyvaerts fc Norman 1989[ 2003. iLivio 20021 



IHonda fc Honda 20021 |Tsinganos fc Bogovalov 2002| ). In 
the relativistic limit however it is slower due to the decolli- 
mating effect of the electric force and the higher inertia of 
the flow, but still possible (Vlahakis & Konigl 2003, 2004). 
Alternatively, collimation in relativistic jets may be due 
to the external pressure of a surrounding slower and easily 



collimated disk wind, Bogovalov & Tsinganos (2005) 



Magnetized non relativistic jets from extended accre- 
tion disks were first modelled analytically by Blandford 
& Payne (1982), wherein the plasma acceleration relies 
on the magnetic extraction of angular momentum and ro- 
tational energy from the underlying cold Keplerian disk. 
This energy is channeled along the large-scale open mag- 
netic ficldlines anchored in the corona or the rotating disk. 
The ionized fluid is forced to follow the fieldlines and to 
rotate with them while it is magnetocentrifugally acceler- 
ated if the angle between the poloidal magnetic field and 
the disk is less than 60°. Cao fc Spruit (1994)| showed that 
in the relativistic case this condition is less severe and that 
close to the black hole the magnetocentrifugal acceleration 
may be efficient even at higher angles. Analytical, radially 
self-similar disk-wind solutions were extended to special 
relativistic cold winds in Li et al. (1992), and Contopoulos 
(1994), by neglecting gravity to allow the separation of the 
variables. Thermal effects were introduced into these rel- 
ativistic models by Vlahakis & Konigl (2003) to analyze 
the formation of a relativistic flow from hot magnetized 
plasmas, showing that such solutions could be applied to 
Gamma Ray Bursts wherein the flow is thermally driven 
at the base. However, most of the acceleration is of mag- 
netic origin and there is an efficient conversion of Poynting 
to kinetic flux of the order of 50%. They also applied this 
disk wind solution to AGN jets ( [Vlahakis fc Konigl 2004| ) 
showing that they could trace the observed parsec scale 
expansion of the wind. Another approach to solve the rel- 
ativistic MHD equations for outflows around black holes 
is to solve numerically the transfield equation in the force- 
free limit IjGamenzind 198fia|l . a study further developed 
in GRMHD by using first a Schwarzschild metric and then 
extending it to a Kerr metric IjFendt 1997(1 . 

Radial self-similarity is usually used in disk- wind mod- 
els due to the complexity of the non linear system of MHD 
equations. However, such solutions cannot describe the 
flow close to the rotational axis where they become singu- 
lar. On the other hand, meridional self-similarity provides 
a better alternative to study the outflow close to the sym- 
metry and rotation axis of the central corona. In the cen- 
tral part where the thermal energy is rather high, the wind 



may be thermally driven. Spherically symmetric relativis- 
tic hydrodynamical models have been proposed to study 
the formation of s uch outflows l(Michel 19721 IDas 19991 
IMeliani et al. 2004)1 . Those models are restricted to the 
case where the magnetic effects in the acceleration are neg- 
ligible. In these models, a wind forms in the hot corona be- 
cause of the internal shock maintained by the centrifugal 
barrier (Ch akrabarti 19891 IDas 200T|l . or by the pressure 
induced via a first order Fermi mechanism IjDas 1999|l . 

An important alternative to analytical models are 
numerical simulations. In the special relativistic do- 
main, simulations have been presented for coronal winds 
(Bogovalov & Tsinganos 1999, Tsinganos & Bogovalov 
2002), and in the general relativistic domain for disk 
winds (e.g. IKoide et al. 19991 2001) to model the for- 
mation and collimation of relativistic outflows in the 
vicinity of black holes. The difficulty for relativistic 
outflows in a single-component model to be collimated 
led Bogovalov & Tsinganos (2005) to propose a two- 



component model wherein a relativistic central wind is 
collimated by a surrounding non relativistic disk-wind. 
Shocks may also develop as the disk-wind collides and 
collimates the inner relativistic wind. Note that all such 
simulations are performed by using time-dependent codes. 
Analytical models conversely have presented more sophis- 
ticated steady solutions of outflows to be used as initial 
conditions in more complex simulations, albeit sacrificing 
freedom on the chosen boundary conditions. 

In this article we present an extension of the 
non relativistic meridionally self-similar solutions 
USauty fc Tsinganos 1994| hereafter ST94) to the case 
of relativistic jets emerging from a spherical corona 
surrounding the central part of a Schwarzschild black 
hole and its inner accretion disk. 

We will not discuss here the origin of the plasma, as- 
suming that it can come from, e.g., the accretion disk or 
pair creation. Furthermore, only the outflow process is 
considered, with the base of the corona placed at a few 
Schwarzschild radii, just above the so called separating 
surface ( Takahash i et al. 19 90 ) . 

Attention is also given to the contribution of the dif- 
ferent mechanisms, hydrodynamic and magnetic, to the 
acceleration and collimation of the outflow, as in previ- 
ous papers of this series dSauty et al. 19991 2002, 2004, 
hereafter STT99, STT02 and STT04). This is also a 
way to extend to 2D outflows, thermally driven, spher- 
ically symmetric (ID) wind models ((Price et al. 2 003 
IMeliani et al.10041). 

In the following two sections, the basic steady axisym- 
metric GRMHD equations are presented using a 3+1 for- 
malism (Sect. 2), together with their integrals. The as- 
sumptions leading to the self-similar model are presented 
in Sect. 3. In Sect. 4, the analytical expressions of the 
model together with the derivation of an extra free in- 
tegral controlling the efficiency of the magnetic rotator, 
as in STT99, are given. In Sect. 5, an asymptotic anal- 
ysis of the solutions is performed as well as the link to 
the boundary conditions in the source. Sect. 6 is devoted 
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to a parametric study of various solutions to emphasize 
the main difference obtained with relativistic flows. We 
discuss the acceleration and collimation of these new so- 
lutions (Sect. 7) and compare them with the non relativis- 
tic model in Sect. 8. In the last Sect. 9 we summarize our 
results and shortly outline their main astrophysical im- 
plications. The confrontation of the present model with 
observed jets from radio loud extragalactic jets, such as 
those associated with FRI and FRII sources, is postponed 
to a following paper, as it involves special techniques for 
constraining the parameters by using observational data 
and a specific iterative scheme to use the model. 

2. Basic equations 

In this section we briefly present, in order to establish 
notation, the governing equations for magnetized fluids in 
the background spacetime of a Schwarzschild black hole 
and also the corresponding MHD integrals. 

2.1. The 3+1 formalism for steady flows 
2.1.1. Schwarzschild metric 

The gravitational potential due to the matter outside the 
black hole is assumed to be negligible. In Schwarzschild 
coordinates (ci, r, 8, ip) the background metric is written 

as, 



ds 2 



where 



-h 2 c 2 dt 2 + Adr 2 



1 

h 2 ' 



2 d 2 9 + r 2 sin 2 <9d^ 2 



1 - 



2QM. 



1 - 



(1) 



(2) 



is the redshift factor induced by gravity at a distance r 
from the central black hole of mass A4», expressed in terms 
of the Schwarzschild radius tq — 2QAi,/c 2 . Note that the 
time line element or the lapse function is usually denoted 
by ho or a. In this paper, for further convenience and 
simplification we have used the symbol h. 

In the following we find convenient to use a 
3+1 split of space-time, following the usual ap- 
proach of MHD flow treatment in general relativ- 
ity IjThorne fc McDonald 19821 IThorne et al. 19861 
|Mobarry fc Lovelace 19861 ICamenzind 1986a|> . The 3+1 
approach allows to obtain equations similar to the 
familiar classical equations. We write all quantities in the 
FIDucial Observer frame of reference, known as FIDO, 
which corresponds to observers in free fall around the 
Schwarzschild black hole. For the FIDO, space time is 
locally flat. 



2.1.2. Particle conservation 

The equation of conservation of particles (nu a )- a = in 
the 3+1 formalism is 

V- (hjuV) =0. (3) 



Here n is the proper number density of particles, V is the 
fluid three-velocity as measured by FIDOs, u a = (jc, r yV) 
and 



7 = (i - v 2 / c 2 r i/2 

is the Lorentz factor. 



(4) 



2.1.3. Maxwell's equations and Ohm's law 

Maxwell's equations written in the 3+1 formulation 
(Thorne & MacDonald 1982, Breitmoser & Camenzind 
2000) are 1 



V • E = 47rp e , 
V-B = 0, 
V x (HE) = 0, 

V x (hB) = — J e , 



(5) 
(G) 
(7) 

(8) 



where (E,B) is the electromagnetic field and (p e , J e ) the 
associated charge and current densities. Ohm's law for a 
plasma of high conductivity is 



E 



V x B 



(9) 



2.1.4. Euler Equation 



Euler's equation is obtained by projecting the conser- 
vation of the energy-momentum tensor, T a ^ = 0, onto 
the spatial coordinates (a = 1,2,3) and combined with 
Maxwell's equations (see IBreitmo ser fc Camenzind~2 000 
or, for an expression closer to ours, albeit restricted to spe- 
cial relativity Goldreich & Julian 1970, Appl & Camenzind 
1993, Heyvaerts & Norman 2003), 



jn(V ■ V) 



(jwV\ = _ Vlnh _ vp 



-1 

-peE 



J P x B 



(10) 



where P and w are the pressure and enthalpy per particle, 
respectively. The form of Eq. 110(1 is equivalent to that in 
Mobarry & Lovelace (1986). 

2.1.5. Thermodynamics 

The first law of thermodynamics is obtained by project- 
ing the conservation of the energy-momentum tensor along 
the fluid four-velocity, u a T^ = 0. In fact, for ideal MHD 
fluid the corresponding contribution of the electromag- 
netic field is null due to the assumed infinite conductivity. 



1 Apparently, Eqs. 15I8H have some differences with their 
equivalent forms in a globally flat spacetime. The h appearing 
in Eqs. J7J and (|SJ results from formulating the laws of differ- 
entiation on a "per-unit-t-base" while the vectors E and B are 
defined on a "per-unit-FIDO-time" IThorne et al. 198611 . 
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Thus, only the thermal energy affects the variation of the 
proper enthalpy of the fluid, 



nV ■ Vw = V ■ VP . 



2.2. Integrals of axisym metric MHD outflows 



(11) 



Assuming axisymmetry of the plasma flow allows us to 
reduce the number of differential equations by integrat- 
ing some of them and thus obtaining conserved quan- 
tities along the streamlines. We follow the notations of 
Tsinganos (1982). 

From Eqs. (JSJ) and © we can introduce a magnetic 
flux function A, 



Bp = V x 



.4 



x e 



(12) 



and a stream function W which gives the particle flux, 



Anh^nVp = V x 



r sin ti l w 



(13) 



V ^ A B 

p 4nhjn p 

Since B p ■ V$/h 



Bp ■ E 



(14) 



B ■ E = 0, the surfaces 



of constant electric potential are also surfaces of constant 
magnetic flux, so $ = &(A). Thus fi = — cd$/dA is also 
a function of A. From Eq. JJjJ the toroidal components V v 
and B v are related 



= z — ; — B<p 



~h~ 



(15) 



This is called the isorotation law because each 
stream/fieldline rotates rigidly with an angular speed f2 
corresponding to the angular speed of the footpoints of 
this poloidal stream/fieldline. 

The azimuthal component of the momentum equation 
yields the conservation of the total specific angular mo- 
mentum, 



w zoB v 
c 2 v V A 



(16) 



The generalized Bernoulli integral (including rest 
mass) 



zuCl 

£ = hjw — h— — B v , 



(17) 



may be obtained by integrating the equation of motion 
along each streamline. The first part of the r.h.s. repre- 
sents the hydrodynamical energy flux transported by the 
fluid while the second part corresponds to the Poynting 
flux. 

We have obtained the usual four integrals of mo- 
tion ^,4, O, £, L ( |Heyvaerts &: Norman 2003| ) that are 
constant along a fieldline for a stationary and axisym- 
metric plasma. They can be used to find algebraic re- 
lations between the Lorentz factor, the toroidal veloc- 
ity and the toroidal magnetic field. Defining the poloidal 
"Alfvenic" number M UMichel 19691 ICamenzind 1986bl 
IBreitmoser <k Camenzind~2 000 ) 



M 2 = 



4Trh 2 nwy 2 V 2 w^ 2 A 



B 2 c 2 



Aimc 2 



(18) 



and the cylindrical distance in units of the light cylinder 
distance (although in our case it is not a cylinder), 

znil 
ch ' 



(19) 



we find, 



where the subscript p denotes the poloidal components 
and w = r sin 9. 

From Eq. JjJ, we can define an electric potential asso- 
ciated to the electric field, E — (V$)//i. Thus, the previ- 
ous equation and axisymmetry imply E v = 0. In addition B^ 
from the flux freezing condition (Eq.|5J) and Eqs. (|12fl and 
Ijl3(l we get that \1/ is constant on surfaces of constant A 
on which the corresponding streamlines and fieldlines are 
roped, V p || B p . It follows that d^/dA = is a function 
of A and we can write 



M 2 x\ 



r: 2 h 2 (1 



M 2 -h 2 (l- x\) 



A 



hjw = £ 



cx 
M 2 



M 2 -h 2 {l-x 2 ) 
h 2 {l-x\ 



M 2 



h 2 (l-x 2 ) 



(20) 



(21) 



(22) 



(23) 



The quantity x\, 
2 _ flL 

is a measure for the amount of energy carried by the 
electromagnetic field. It is a measure of the energy flux 
of the magnetic rotator in units of the total energy flux 
( Breitmos er fc Camenzindl2 000 ) . 

The MHD equations possess a well known singularity 
at the Alfven surface where the denominator of Eqs. (|20() . 
(|21|l and (|22|l vanishes. Then, the numerators should van- 
ish simultaneously to ensure a regular behavior, implying 
at the Alfven point, 

(24) 



I Alfven 



h*c 



X A ■ 



M 



Alfven 



= K l-xl) =Mi 



Using Eq. ijT%|) we find 



B 



2^2 



47T7^n*ui + 



(25) 



(26) 



where the subscript 7k- denotes quantities evaluated at the 
Alfven point along the polar axis. Note that the position 
of the Alfven surface is shifted with respect to the classical 
case because of the lapse function h and the existence of 
the light cylinder x = 1. 

In the Newtonian limit, Eqs. J2U and give L = 
mvjpQ. where m = w/c 2 is the particle mass. 
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3. Construction of the meridionally self-similar 
model 

Our goal in this section is to find semi-analytical solutions 
of the r— and 9— components of the Euler Eq. H10|) . by 
means of separating the variables r and 9. In order to 
facilitate the analysis it is convenient to use dimcnsionlcss 
quantities normalizing at the Alfven radius along the polar 
axis. Using the notations introduced in ST94 we define a 
dimensionless radius R and magnetic flux function a, 



R 



r 



.4: 



-a . 



(27) 



We introduce two dimensionless parameters to de- 
scribe the gravitational potential. The first is v which 
represents the escape speed at the Alfven point along the 
polar axis in units of V*, 



I2GM. 



(28) 



The second parameter 2 is the ratio of the Schwarzschild 
radius over the Alfven radius r*. 



Vi 



(29) 



which is also the escape speed in units of the speed of 
light. 

Combining Eqs. I|28l) and (|29|l we get a condition that 
restrict the parametric space to 



(30) 



3.1. Separation of the variables 

3.1.1. Magnetic flux and "Alfvenic" number 

As a = on the rotational axis and we are interested 
on the central component of the jet, we assume that the 
cross section area of a given magnetic flux tube can be 
expanded to first order in a, 

S(R, a) = tttd 2 = irrlG 2 (R)a . (31) 

Normalizing at the Alfven surface, we choose G(R = 1) = 
1 such that G is the cylindrical radius in units of the 
Alfvenic cylindrical radius. Thus G(R) = w/uj^ with 

This is equivalent to assume, as in ST94, that the mag- 
netic flux function a(R, 0) has a dipolar latitudinal depen- 
dence, 



R 2 



■ sin 2 9 . 



G 2 {R) 

We also introduce the expansion factor F 



F = 



dim 



d\nR 



dlnG 2 
dhiR 



(32) 



(33) 



2 u corresponds to the parameter m used by Daigne & 
Drenkhahn (2002). 



In addition, we assume that the surfaces of constant 
"Alfvenic" number are spheres, such that, 



(34) 



M 2 (R,a) = M 2 (R) , 
with M 2 (R = 1) = hi. 

3.1.2. Pressure 



The pressure dependence is obtained by making a first 
order expansion in a 



I 7 ,V^v; 2 n(i?)(i + Ka) 



r = P a + -^n^V^U(R)(l + K a), (35) 
with P , k constants and II (R) a dimensionless function. 



3.1.3. Free integrals 

Combining Eqs. (TTH|l and l(3"4T we deduce that 
47rnc 2 M 2 



Expanding the r.h.s. to first order in a, we find 



47rc 2l^± (1+fo) 



(36) 



(37) 



where 8 is a free parameter describing the deviations from 
spherical symmetry of the ratio number density/enthalpy 
and not of the density itself as in ST94. 
Similarly we expand a 



(38) 



where A is a constant measuring rotation. 

Finally we choose for a form similar to the one in 
ST94 



Vl + 5a 



(39) 



From Eq. I|19|l we can now express x 2 in terms of a. 
and G{R), 



V 2 h 2 
\ 2 \^±G 2 



c 2 h 2 1 + 5a 
Similarly, from Eq. 124|) we find 



M^ 2 hi 2 
v 2 h 2 1 



5a 



T 2 

x A 



V 2 



jiX- 



(40) 



(41) 



which gives from Eq. Il23[l the form of the Bernoulli inte- 
gral, 



(42) 



Note that the parameter A in Eq. (|39|l is the same constant 
as in Eq. (|3*5|l because the energy £ must be equal to its 
hydrodynamic part h^^w^ on the rotational axis where 
the Poynting flux vanishes. 
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3.2. Electric force 

Conversely to the non relativistic limit, we cannot neglect 
the charge separation and the presence of the electric field. 
From the previous assumptions, we can calculate the elec- 
tric force p e E which has the following two components 

r2 ( l2 jp 

' ' : x\G 2 sm 2 6 



p e E r 



Airr+G 4 1 h 2 



4 



'dR 



El 

R 



2F- 



R(l + Sa) 

3±L x 2 %r^ 2 9^ 

c 2 (1 



h 2 R 



8a 



PeE 



x\ G 2 sin 9 cos ( 




R(l+Sa) 

2 V 2 2 
A 2 -4rR sin 9 cos 9 " 



8a 



(1 + Sa) 2 



(44) 



Note that in this form of the expressions of the forces the 
variables are not separable. 

3.3. Non relativistic rotation 

Contrary to the non relativistic case and in order to sep- 
arate the variables (R, 9) in the r— and 6— components 
of the Euler Eq. l|l(Jf) we need some further assumptions. 
Basically we expand these equations with respect to 9. 

We suppose that the rotational speed of the fluid re- 
mains always subrelativistic. In other words, we focus on 
streamlines that never cross the light cylinder such that 
the later does not affect the dynamics (x -C 1, which im- 
plies a; a <C 1). Of course this reduces the domain of valid- 
ity of the solutions to the vicinity of the rotational axis be- 
cause x\ as it is given by Eq. I|41|) is sufficiently small only 
for relatively small a. The region of validity of our model 
depends on how small the parameter X 2 V 2 /c 2 = X 2 p,/v 2 
is, though. The requirement that the light cylinder lies fur- 
ther away from the Alfvenic surface (xa < 1) constrains 
the parametric space to S > \ 2 p/v 2 . On the other hand, 
there is a possibility that all streamlines never cross the 
light cylinder. This happens when the equation 



1 O 



1 _ AV G 2 hj 
a v 2 h 2 



(45) 



cannot be satisfied for any R, or equivalently 8 > 
{\ 2 p/v 2 ) [G 2 h 2 /h 2 '\ . Weak rotation (small A) or sig- 
nificant deviation of the particle flux from spherical sym- 
metry (high 8 which results in a fast decrease of fl as we 
move away from the rotation axis) contribute to the va- 
lidity of the above inequality. 



A consequence of this is that the jet is thermally 
driven. Indeed the ratio between the Poynting flux and 
the matter energy flux is, 



-hmVLB v /-fy A 



= h 2 



(46) 



(43) B, 



h-fw " M 2 - h 2 {\- x\) 

Thus, in the approximation x -C 1 the contribution of the 
Poynting flux is negligible in accelerating the flow. 

More generally, after expanding with respect to x A , we 
neglect terms of the order x\sva.9 or higher. This is also 
consistent with our previous assumption of keeping terms 
only up to a in the integrals. For example, the expression 
of the azimuthal magnetic field, Eq. (|21|l . becomes 

\B± K G 2 hl - h 2 



-Rsin9. 



(47) 



G 2 h M 2 -h 2 

Equivalently, after expanding all terms with respect to 9 
we neglect terms of the order of sin 3 9 or higher. 

As another example, the approximate form of the elec- 
tric force is 

^«7Z^^A 2 f W0, (48) 



B 



Airr+G 4 h 

2 \ 2 V 2 

— A —7T R sin 



P ° E °* 2 4^Gl{T) ^i^cos*, (49) 

where we have further approximated T^jj^p ~ 2, since 

this factor is multiplied with x\/ sm9 cx sin 9 and a cx 
sin 2 9. Though this approximation is consistent with the 
fact that we neglect any terms of the order of sin 3 9 or 
higher, it gives an extra restriction. Thus, the model ap- 
plies only to the region near the rotational axis where 
a <g; 1/8. We shall calculate both inequalities (x A <C 1 
and a <C 1 / 8) a posteriori in order to determine the regime 
where the solution is valid. 

Note that in the very vicinity of the black hole h — > 
and since the factor hx does not vanish (see Eq. I19fl , x be- 
comes larger than unity, as expected by the presence of the 
second light cylinder close to the horizon (e.g., Takahashi 
et al. 1990). However, it is enough that the coronal base is 
at a distance of a few gravitational radii (R>2p) to guar- 
antee that x <C 1 at the base of the outflow (this condition 
is fulfilled in our solutions). 

4. Equations of the model 

4.1. Expressions of the fields and the enthalpy 

The velocity and magnetic fields can now be written ex- 
clusively in terms of unknown functions of R. For later 
convenience, as in ST94, we shall denote by Nb, Ny and 
D the following quantities that appear in various compo- 
nents of the fields, 

N B = S - G 2 , (50) 



AV = 



D 



G 



h 2 
M 2 

/i 2 _ M 2 
h 2 hi 



(51) 
(52) 
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Thus, 



cos 9 . 



sin ( 



r G 2 

B* hF 

Be CP~2 

XB* K N B 
* ~ G 2 h D 



R sin 9 , 



V r 



KM 2 



1 



h 2 G 2 VTTfe 
KM 2 hF 1 



/jA 2 7V £ 
j/ 2 D 



•Q 



/i 2 g 2 2 yrr^ 

/i AK AV i?sin6> 



sm ( 



V * KG 2 D y/i 



5a 



(53) 
(54) 
(55) 
(56) 
(57) 
(58) 



The electric field can be deduced from the flux freezing 
condition and the above equations (see Eqs. BBlBSt - 

Similarly, the enthalpy and the particle number density 
are given by 



hjw = h+~f*w+ 1 



v 2 D 



-a 



hjn = a*7*7i*-^ [ 1 + da — -ppa 



M 2 



v 2 D 



(59) 



(60) 



where we used Eq. IpTfjl to deduce Eq. ijfinjl , while the pres- 
sure is given by Eq. |3"B*)l. 

4.2. Ordinary differential equations and numerical 
technique 

There are three equations given in Appendix [S] that de- 
termine the three unknown functions H(R), F(R) and 
M 2 (R). We recall that G 2 is related to F through Eq. 

Before discussing in detail the results of the para- 
metric study we outline the method for the numerical in- 
tegration of Eq. (JsSJ and Eqs. (jA.ljl - (|A.3jl . We start 
integrating the equations from the Alfven critical surface. 
In order to calculate the toroidal components of the fields, 
i.e. N B /D and N v /D = N B /D - 1, we apply L'Hopital's 
rule, 



Nb 
D 



hl(2 - F+) - fi 
hip 



1 dM 2 



h 2 dR 



-a- (61) 



To avoid kinks in the fieldline shape, we need to satisfy 
a regularity condition ( |Heyvaerts fc Norman 1989| ). This 
means that Eq. I|A.3|) should be regular at R = 1. As in 
ST94 this extra requirement is equivalent to Afp,D = 
which eventually gives a third order polynomial equation 
for p 



C 3P 3 + C 2 p 2 + C lP + C o = 0.. 



Co = -— F* - 2 



h 2 



C 1 = A^-2 + t ], 



(62) 

(63) 
(64) 



C 2 

c 3 



= hl^ 



^F 2 + i(i- K n.) + A 2 4 



(65) 
(66) 



Once we have determined the regularity conditions at 
the Alfven point, we integrate downwind and upwind and 
cross all the other existing critical points as in the non 
relativistic case. 

It is worth noticing that the shape of the streamlines 
F± at R = 1 is determined by the regular crossing of the 
slow magnetosonic surface. We point out further that, be- 
sides the free parameters listed at the beginning of the 
section, solutions depend also on IT*, i.e. the pressure at 
the Alfven surface. As in the classical case its value has 
been chosen such that the gas pressure is always positive. 
More details on the numerical technique can be found in 
ST94 and STT02. 

4.3. The integral e 

As in ST94, it is possible to find a constant e for all field- 
lines. This parameter e has been used in ST94 and in the 
following papers to classify the various solutions. We shall 
use a similar technique to construct such a constant in the 
present model. 

Eq. after substituting n from Eq. (|18|l and using 
V ■ V cx d/dR\ a (derivative keeping a constant), can be 
re-written as 



8ttM 2 



dP 
dR 



d_ 

~dR 
d_ 
dR 



^ 2 A (E 2 -w 2 )^ 



(67) 



where ^> 2 A w 2 is proportional to the energy per unit volume 
of the fluid in the comoving frame, i.e. reduced to the 
thermal content. Thus ^\(£ 2 — w 2 ) in essence measures 
the variation between the total energy and the thermal 
energy of the fluid. 
By writing the term 



^ 2 A w 2 



*> 2 7 2 (1 



V 2 /c 2 -V 2 /c 2 ) 



M 4 B 2 
~h 2 ~ 



(68) 



and using Eqs. lf2TJ|) and J22J) we find 



^ 2 A w 2 



V 2 A £ 2 
c 2 h 2 



M 2 -h 2 (l-x 2 A ) 
M 2 -h 2 {l-x 2 ) 



V 2 A £ 2 



2 ? 2 tM 2 x 2 Jx-xh 2 {\-x 2 A ) 



c 2 h 2 



M 2 -h 2 {l-x 2 



M 4 B 2 
~H 2 ~ 



(69) 



In the particular model that we examine, the form of 
the pressure is P = /i(i?)(l + Ka)/8-7r. We also know 
the 9 dependence in all quantities in the expression for 
w 2 /c 2 , and after expanding with respect to sin 2 9 we 
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find ^ 2 A (£ 2 
gives 



w 2 )/c 2 = f 2 {R) + h(R)a. Then Eq. J2) 



M' 



:dfl 

dR 



dR 



df 3 



dR 

-M 2 dh = d/ 2 
—M 2 ndfi = d/ 3 



(70) 



Eliminating dfi we get the integral — nfziR) = 

e =const. After substituting the expressions for /2CR) and 
fa(R), we arrive at 



Af 4 



hfR 2 G 2 



El 
4 



1 

ft 2 



I? 2 



/i 2 G 2 
D 



(5 - k) v 2 
h 2 R 

2A 2 N B 



where e is a constant, the same for all ficldlines. 

Similarly to what was done in ST99, we can calcu- 
late this constant at the base of the flow R assuming the 
poloidal velocity is negligible there [M(R ) « 0]. Let's ex- 
press e/2A 2 in terms of the conditions at the source bound- 
ary, 



6 



R,o 



^Poynt.,o 



2A 2 



(72) 



where f mr = h 2 LQ is the energy of the magnetic rotator, 
£poynt. = —hvjQ,B v /ty a is the Poynting flux and 



tp,o 



(73) 



is the rotational energy per particle. It is proportional to 
the specific rotational energy V 2 G /2, with the factor £/c 2 
having the dimensions of a mass. Finally we have 



AS* G = - 



£ /ic 2 (5 — k) a 
~^?2R l+Ja 



(74) 



a term similar to the nonrelativistic case where it mea- 
sures the excess or the deficit on a non polar streamline, 
compared to the polar one, of the gravitational energy per 
unit mass which is not compensated by the thermal driv- 
ing (STT99). As in the classical case, e measures the effi- 
ciency of the magnetic rotator to collimate the flow. Thus 
if e > we have an Efficient Magnetic Rotator (EMR) 
where magnetic collimation may dominate, while if e < 
we have an Inefficient Magnetic Rotator (IMR) where col- 
limation cannot be but of thermal origin. 

5. Asymptotic behaviour 

In the region far from the base where the jet attains its 
asymptotic velocity, assuming it becomes cylindrical, the 
forces in the radial direction become negligible, since the 
jet is no longer accelerated. In the transverse direction, 
the following four forces balance each other: the transverse 
pressure gradient, fp, the total magnetic stress (magnetic 
pinching plus magnetic pressure gradient) of the toroidal 



magnetic field component, fs, the centrifugal force, fc 
and the charge separation electric force, f e, 



fc + f B + fp + ft 



0, 



(75) 



The full expressions of these forces are given in Appendix 
IbI In the asymptotic region, 9 ~ 0, they can be written 
as follows in cylindrical coordinates, 



fc = 



(71) h 



fp 
f E 




(76) 



2irG! 



ZU 



(77) 
(78) 
(79) 



The centrifugal and electric forces have a decollimating 
effect on the jet, while the pinching magnetic force colli- 
mates it because of our choice on the current distribution. 
For asymptotically underpressured jets where k > and 
TToo > (or k < and Iloo < 0), the pressure increases 
away from the polar axis which helps collimation. The 
opposite holds for overpressured jets where k < and 
Iloo > (or k > and LToo < 0). 

Combining Eq. (^U with Eqs. (JZEJ) - JZHJ we obtain, 



2A 2 



IT 



1 



2M 2 



N 



Voc 



fi 



N, 



Boo 



D, 



(80) 



The second equation controlling the asymptotic trans- 
verse force balance is given by e, Eq. 1)71(1. In the asymp- 
totic region, for a cylindrically collimated jet — > 2 and 
e becomes, 



e 
2A 2 



Mi 



1 



(81) 



We notice that the asymptotic behaviour of the jet is de- 
scribed by the asymptotic pressure Iloo and the two pa- 
rameters e/2A 2 and k/2X 2 . These equations are similar to 
the classical model except for the decollimating effect of 
the electric field and charge separation. Besides that note 
also that in D, Ny and Nb the space curvature at the 
Alfven critical surface also appears. 

At the base of the wind the Alfven number vanishes, 
M — » 0, while the opening of the jet is weak, G 2 <§C 1. We 
can use this criterion to define the distance R a where the 
outflow starts. Namely from the expression of e, Eq. (|71|) . 
we get, 



Ro 



(<5 — k)v 2 — fie 
2A 2 -e 



(82) 
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We see that in order to have acceleration (R a > fj,) we 
must have approximately (S— k)/2 > /.iA 2 jv 2 for 2A 2 — e > 
0, extending the criterion found in the classical regime 
(8 > k; the above expression reduces to Eq. 14 of STT02 
for \i — > 0). In particular we note that for e > or e < 
relativistic effects enlarge or reduce the size of the sub- 
Alfvenic region, respectively. 

Note that the definition of R a coincide with the so 
called separating surface (see e.sr. lTakahashi et al. 199fll for 
a cold plasma in a Kerr metric). Above this surface the 
plasma is outflowing. Below this surface other critical sur- 
faces exist (e.g., IBeskin fe Kuznetsova 2000jl but remain 
out of our consideration. If the plasma is created via pair 
production this may constrain the boundary conditions at 
the base of the flow. However we do not discuss the origin 
of the coronal plasma in this paper. 

6. Parametric analysis 

As as first step of the numerical analysis we have per- 
formed a study of the effects on the solution of the free 
parameters of the model /i, v, A, S and n. With respect to 
the classical case, the relativistic effects are ruled by the 
new parameter fi. 

6.1. Effect of space curvature and gravity (fi, v) 

The parameters \i and v denote the escape speed in units 
of the light and Alfven speed, respectively. However simi- 
lar they look, they have opposite effects on the initial ac- 
celeration and the terminal speed. In the super- Alfvenic 
region the acceleration is not strongly affected by differ- 
ent values of /i for /i < 0.1; in fact, the effect of relativistic 
gravity is negligible after lOrg (Fig. QJ. So the effects we 
are discussing now refer to the base of the flow, in the 
sub Alfvenic regime. 

The parameter /i, the ratio between the Schwarzschild 
and the Alfvenic radius, representing also the escape speed 
at the Alfven radius in units of c, is related to the relativis- 
tic effects of gravity in this model. Basically fi controls the 
extension of the corona and the acceleration of the flow 
in the sub- Alfvenic region. Increasing fi increases also the 
asymptotic velocity, as it can be seen in Fig. A simple 
physical interpretation may be given to this behaviour. 
When the Alfvenic surface approaches the Schwarzschild 
surface, gravity in the sub-Alfvenic region, and thus in 
the corona, increases. Consequently, to support gravity 
the thermal energy increases too. This larger amount of 
thermal energy in the corona will be transformed in turn 
largely into kinetic energy along the flow. In other words, 
the increase of y, implies a stronger density gradient of the 
flow in the sub-Alfvenic region, increasing the radial pres- 
sure gradient dH/dR and leading to a stronger expansion 
and acceleration. 

This behaviour can be also understood considering 
that larger values of \i imply a larger space curvature, 
increasing also the expansion of the streamlines and 
thence the efficiency of the acceleration, as it has been 



shown in the study of the relativistic Parker wind (see 
IMeliani et al. 2004jl . 

Conversely, increasing v decreases the asymptotic ve- 
locity as well, since it reduces the size of the corona, keep- 
ing fi constant, that is the distance of the Alfven surface 
to the Black Hole (see Fig. ^3). This figure shows that 
the base of the flow gets closer to the Alfven radius and 
farther from the Schwarzschild radius. As in the non rela- 
tivistic case the parameter v 2 is the ratio of gravitational 
to kinetic energy at the Alfven surface, Eq. |2H1 An in- 
crease of v reduces the fluid velocity at the Alfven ra- 
dius with respect to escape speed needed to get out of the 
black hole's attraction. The reduction of the size of the 
corona is also consistent with the reduction of the veloc- 
ity. The behaviour is as expected from the solutions in the 
classical regime (ST94, STT02, STT04): the higher is the 
value of v the lower is the asymptotic velocity, although we 
didn't show it explicitly, as in the present Fig. There 
is also a minimum value of v to have mass ejection, below 
that value the thermal energy cannot support gravity (see 
STT02). 

6.2. Effect of rotation (X) 

The parameter A is related to the rotation of the flow and 
to the axial electric current (Fig. Eb)- As for non rela- 
tivistic outflows (ST94, STT02) it rules the jet dynamics 
through the Lorentz force, collimating asymptotically the 
jet via the toroidal magnetic field, while the centrifugal 
force has instead a decollimating effect. We have checked 
that the behaviour of the solutions is similar to the clas- 
sical case. The increase of A leads to more collimated and 
slower jets (Fig. Et) • This can be understood as follows. 
Increasing A increases the axial current (Fig. |5Jd) which 
increases the toroidal magnetic pinching. In order to pre- 
serve equilibrium the flows reacts by increasing the cen- 
trifugal force and thus its rotational speed by reducing 
its cross section. The reduction of the expansion factor 
reduces as usual the pressure gradient and the thermal 
driving efficiency thus reducing the asymptotic speed. 

In addition we must take into account that A is re- 
lated to the electric potential $. Consequently it controls 
the charge separation and the corresponding electric force 
(p e E oc x 2 oc A 2 ). This force becomes dominant where the 
jet rotation speed becomes relativistic {wVL ~ c) . In other 
words the higher is A the larger the effects of the light 
cylinder (Fig. |3J) . 

6.3. Effect of pressure and density anisotropies (k, 6) 

The physical meaning of k remains the same as in non rel- 
ativistic flows. For k positive or negative the gas pressure 
increases or decreases with colatitude, respectively. Then 
in the first case the gas contributes to the thermal confine- 
ment of the flow (underpressured jets), while in the second 
to its thermal support (overpressured jets). We have lim- 
ited ourselves to present here the study of underpressured 
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Fig. 1. Variation of the velocity vs. r /r* for different values of fi in (a) and v in (b). The other parameters are 
v = 1.185, 8 = 2.113, k = 0.5, A = 1.9995 in (a), and fi = 0.1, 6 = 1.0, k = 0.5, A = 1.0 in (b). 
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Fig. 2. Comparison between two jet solutions with A = 1.0 
(e = 0.089) and 1.2 (e = 0.824). In (a) we plot the velocity 
along the polar axis and in (b) the dimcnsionless electric 
current density along the polar axis. The other parameters 
are ^ = 0.1, v = 0.65, 8 = 1.4, k = 0.2 and IT* = 0.75. 

flows. The behaviour of the relativistic solutions with n is 
analogous to that of classical flows. For higher k both the 
asymptotic velocity and jet cross sections decrease (see 
STT99 and STT02 for details). 

In the present model the parameter 8 controls the vari- 
ation of the ratio n/w with colatitude, or equivalently in 
the direction perpendicular to the rotational axis. This is 
a relativistic generalization of the classical solutions where 
it governs the transverse profile of the mass density. As a 
result, the effects are similar to those found for non rela- 
tivistic solutions. A larger 8 means a larger gravitational 



potential of the external streamlines with respect to the 
axis, where the acceleration is more efficient. Then the 
asymptotic velocity increases with 8 (ST94, STT02). 

7. Jet dynamics, acceleration and collimation 

We will address now the question of the process of accel- 
eration and collimation of the jet in the case of EMRs and 
IMRs. Keeping fixed /i — 0.1 we have displayed the results 
for an IMR solution with v = 0.781, 8 = 2.613, n = 0.490, 
A = 0.880 and IL. = 1.4 (e = -1.747) in Fig. Hi (see also 
Figs.Et andgt), and an EMR with v = 0.541, 8 = 3.253, 
K = 0.39, A = 1.401 and IL. = 1.26 (e = 1.128) in Fig.g^ 
(see also Figs.Et andEb)- 

We see that the shape of the streamlines in the two 
cases (Figs. 0J) looks similar to the corresponding non rel- 
ativistic case (STT02, see later Figs.[7J|. IMRs show a fast 
expansion in an intermediate region, while far from the 
base the collimated streamlines show strong oscillations. 
EMRs show conversely a continuous expansion with rela- 
tively mild oscillations, or even no oscillations at all when 
pressures are lower. We also display for the two solutions 
the energies along a given streamline (Fig. [SJ and the 
forces perpendicular to the flow (Fig.EJ). 

7.1. Acceleration 

By construction of this model, the wind is basically ther- 
mally driven. At the lowest order we have £ ~ h^fw, while 
the first order term corresponds to the Poynting flux which 
remains however of the same order than the thermal terms 
in the transverse direction. This supposes that there is a 
high temperature corona around the black hole as pro- 
posed by Chakrabarti (1989) and Das (1999, 2000). The 
latest has shown that the stronger is the thermal energy, 
the more stable is the corona. 

We can study the acceleration of the jet analysing the 
contribution of the different energies and their conversion 
from one to the other along the streamlines. The dominat- 
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Fig. 3. Morphology in the poloidal plane of the streamlines of the two solutions of the previous figure and their light 
cylinders for A = 1.0 (e = 0.089) in (a) and A = 1.2 (e = 0.824) in (b). As in Fig. [2 the other parameters are fi = 0.1, 
v = 0.65, 8 = 1.4, k = 0.2 and II* = 0.75. The light cylinders are represented by the two thick solid lines which 
surround the jets, and the different regions of validity of our solutions are shown: solid, dashed and dashed-dotted 



lines correspond to streamlines where the two quantities (defined in Sec. 2.2) x\G 2 and (2 
< 0.01, < 0.1 and > 0.1, respectively. 



8a) {I + 8a) - 2 are 




Fig. 4. Morphology of the poloidal streamlines for two solutions corresponding to an Inefficient Magnetic Rotator 
(IMR, e = -1.747) in (a), and an Efficient Magnetic Rotator (EMR, e = 1.128) in (b). We chose \i = 0.1 in both cases, 
while the other parameters are v = 0.781, k = 0.49, 8 = 2.613, A = 0.880 and II* = 1.4 in (a), and v — 0.541, k = 0.39 
, 8 = 3.253, A = 1.401 and II* = 1.26 in (b). The various regions of validity of our solutions are shown as in Fig. |3] 



ing energy at the base of the outflow is the enthalpy. Part 
of it is used to balance gravity and the remaining part is 
converted into kinetic energy in the region of expansion 
of the jet and stops when the streamlines collimate (com- 
pare Figs. 0] and EJ). In fact, during the expansion of the 
jet, the plasma density decreases which also induces a de- 
crease of the enthalpy. In turn, it creates a strong pressure 



gradient (VP = nWw) that accelerates the jet. For the pa- 
rameters we have chosen, we see that the EMR solution 
has a larger expansion factor than the IMR one (Figs. QJ, 
because of the thermal driving. This is correlated to the 
larger Lorentz factor of the EMR solution (7 ~ 2.8) as 
compared to the IMR one (7 ~ 2.4). This result is not 
related to the nature of the magnetic rotator. For other 
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Fig. 5. Variation of the energy flux normalized to the mass energy, along the external streamline for the IMR solution 
in (a) and the EMR solution in (b) of the previous figure. The Poynting energy is < 1% of the enthalpy at the base 
of the flow and is not plotted. The parameters are the same as in Fig. 0] 
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Fig. 6. Plot of the transverse forces for the relativistic IMR (a) and EMR (b). Forces are normalized by their maximum 
value, usually reached at the base of the flow. The parameters are the same as in Fig. 0] 



parameters, we would get different results but the asymp- 
totic Lorentz factor always increases with the increase of 
the expansion factor because of the thermal driving. We 
verified that the Poynting flux in the two solutions is neg- 
ligible, representing at maximum only 1% of the enthalpy 
at the basis of the flow. 

The IMR solution undergoes a strong expansion in this 
region until a distance of lOOr* and then recollimates and 
consequently decelerates because of the compression. On 
the other hand, the EMR solution collimates already at a 
distance of approximately 50?"* and accelerates all the way 
downwind. In other words, the nature of the collimation 
affects obviously the velocity profile. 



7.2. Collimation 

The collimation of the flow is controlled by different types 
of forces that depend on the morphology of the streamlines 
in the jet. We have plotted for the two solutions the forces 
perpendicular to the streamlines in Fig. [H] Asymptotically 
the centrifugal force is the dominant term which supports 
the wind against either the magnetic confinement in EMR 
or the pressure gradient in IMR. 

The behaviour of the other forces depends on the 
shape of the streamlines (Fig. 0) and they play a rele- 
vant role in the intermediate region before collimation is 
fully achieved. In particular, the stress tensor from the 



poloidal magnetic field and the gravity favor deviation 
from radial expansion while the thermal pressure gradient 
initially tends to maintain the radial expansion. In the re- 
gion of formation of the jet, the strong gravity along the 
polar axis generates a strong pressure gradient. As den- 
sity and pressure increase with colatitude (<5, K > 0), it also 
generates a total force /vj_p + fv± lnfc °c ^ 2 (^ — K ) which 
further out in the jet provides the thermal confinement. 

In an IMR, neither the pressure gradient nor the pinch- 
ing force from the toroidal magnetic field can brake the 
expansion of the flow, and the recollimation occurs where 
the curvature of the poloidal streamlines becomes rele- 
vant. In the asymptotic region, the collimation is mainly 
provided by the transverse pressure gradient, /til, which 
balances the centrifugal force. The pressure confined jets 
undergo strong oscillations similarly to the non relativistic 
case (STT94, STT99). 

In an EMR, conversely, the pinching force of the 
toroidal magnetic field provides collimation all along the 
flow. The pressure gradient may help this collimation as 
for the present solution or be completely negligible for 
other sets of parameters. The magnetic pinching force is 
balanced asymptotically mainly by the centrifugal force 
which tends to decollimate the jet. We must notice also 
that, as expected in the relativistic case, the electric force 
is always positive and its effect in decollimating the jet 
may be comparable with the centrifugal force, differently 
from IMRs (see FigEJ). 
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Fig. 7. Morphology of the poloidal streamlines for non relativistic IMR (a) and EMR (b): the parameters are the 
same as in Fig. 4. with fi = 10~ 5 . 




Fig. 8. Plot of the transverse forces for the non relativistic IMR (a) and EMR (b). We assumed /i — 10 5 while the 
other parameters are identical to the corresponding relativistic solutions displayed in Figs.fJ] 



8. Relativistic vs. non relativistic outflows 

We analyse here in more detail the main differences be- 
tween the relativistic and non relativistic solutions already 
discussed in previous Sect. 6. By increasing the value of \i 
we increase the depth of the potential well. In our calcula- 
tions we have assumed /i = 0.1 for the relativistic solutions 
and this can be justified as follows. 

We supposed that the Alfven surface is roughly at a 
distance of 10r o from the central object. This typical dis- 
tance is usually chosen because it corresponds to the case 
where the wind carries away all the accreted angular mo- 
mentum, provided about 10% of the accreted mass goes 
to the jet ULivio 1999JI . This is of course arbitrary but al- 
lows us to compare our solutions to other models. In the 
case of young stellar jets, the star has a mass of the or- 
der of ~ 1M Q . The corresponding Schwarzschild radius is 
approximately tq w 3 km. Therefore, space curvature at 



the Alfven surface corresponds to a value of ji ss 10~ 5 . 
Conversely, for AGN jets, with a central super massive 
black hole of ~ 10 9 Af Q , we have tq ~ lO 4 i?0 which cor- 
responds to the value we have chosen /i w 0.1. 

We have drawn the corresponding morphologies of two 
non relativistic solutions associated with an IMR and an 
EMR in Figs.0 keeping the other parameters as in Fig. 
0] We also plotted the corresponding transverse forces for 
the non relativistic solutions in Figs. |H1 

8.1. Jets from Inefficient Magnetic Rotators (IMR) 

Let's first turn our attention to the solution from an IMR. 
The morphologies of classical and relativistic jets show in- 
deed small differences. The relativistic jet, though, under- 
goes an expansion slightly more important in the interme- 
diate region (Fig. Et) than in the classical case (Fig. Et). 
This expansion induces a slight relative increase of the cur- 
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EMR 










IMR 




Log(l/V ) 








(b) : 


1 IMR 


EMR 





Log(l/[L ) 

Fig. 9. Plot of the asymptotic dimensionless cross section 
of the jet G 2 ^ in (a), and the asymptotic velocity in (b) 
as a function of the parameter /i on a logarithmic scale. 
The other parameters are those given for the previous rel- 
ativistic IMR and EMR solutions. 



vature forces (inertial and magnetic) compared to other 
forces. In the asymptotic region, the relativistic jet rec- 
ollimation is comparable to the non relativistic one (Fig. 

Note also that in the asymptotic region, the relativistic 
jet pinching by the toroidal magnetic field is almost null, 
while in the non relativistic solution this force is of the 
order of the pressure gradient (Figs. EH and EH)- This be- 
haviour is a consequence of the decrease of the collimation 
efficiency in relativistic jets. 

Last, we know that in the relativistic solutions, there 
is an extra electric force, p e E, due to the non negligible 
charge separation, which also decollimates. In IMRs, its 
influence remains however weak because of the low mag- 
netic field. 

8.2. Jets from Efficient Magnetic Rotators (EMR) 

In the jet solution associated with an EMR, the solution is 
effectively magnetically collimated because of the toroidal 
magnetic field pinching force. This solution clearly shows 
the role of the magnetic field in collimating the jet but also 
the contribution of the charge separation to the electric 
field and then to the decollimation. 

The morphology of EMR jets is more affected by rela- 
tivistic effects than the one of IMR jets. This can be seen 
by comparing Fig. and Fig. [TJd. In particular, the jet 
radius, or equivalently the expansion factor, asymptoti- 
cally becomes more important in relativistic jets because 




Classical 



Relativistic 



r/r 



Fig. 10. Density of the electric current normalized I z , in 
the classical and the relativistic jets for EMR shown in 
the previous figures. 



of the increase of the centrifugal and electric forces (Figs. 
Et> andEt). Simultaneously, the magnitude of the Lorentz 
force decreases and the thermal acceleration increases. 

We can give a simple explanation to this relativistic 
effect. An increase of gravity at the base of the jet in- 
duces a decrease of the normalized electric current in the 
jet because / cx h* where / = I /(c/2r+B+). We have plot- 
ted the dimensionless electric current in Fig. (|10|l for the 
two EMR solutions. The electric current I flows through 
a given cross sectional area S = irm 2 of the classical 
and relativistic solutions. We use the normalized electric 
current because of the different scaling of the classical 
and the relativistic solutions. The decrease of the cur- 
rent goes with a decrease of the toroidal magnetic field 
and with an increase of the expansion of the jet as ex- 
plained in Sec. 16.21 and, consequently, with an increase of 
the poloidal velocity. Note that this increase of the ve- 
locity corresponds to the increase of the relativistic grav- 
ity as we already discussed. The new point is that, con- 
versely to hydrodynamical models, it also decollimates the 
flow. On the other hand, the magnetocentrifugal driving 
of the Poynting flux becomes weaker in relativistic ther- 
mally driven winds, as expected. Thus, as the rest mass 
increases, the Poynting flux is getting weaker relatively to 
the other energies, £p ynt. <C mc 2 , while, the thermal en- 
ergy becomes relativistic and comparable to the rest mass 
(w — mc 2 ) « mc 2 . 

The relativistic effects on the jet acceleration become 
remarkable only for a distance between the Alfven surface 
and the Schwarzschild surface smaller than 100 (ic. /i > 

0. 01, see Fig. HTk,). For fi < 0.01, using this model, we 
see that outflows from a star with mass 1-M© and starting 
at lOOrc are simply scaling with n fFigs. ITTb and II lb') . 
Similarly the ratio between the energetics of the two jets 
are simply proportional to /i. The scaling just reflects the 
linear growth of the flow formation region with gravity, 

1. e. with /x. 

Conversely, for /i > 0.01, the jet is formed at a dis- 
tance smaller than lOOr^, this linear scaling with fj, of 
the dynamics and the energetics does not hold any longer 
because of non linear relativistic effects. The thermal en- 
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ergy converts more efficiently into kinetic energy (Fig. lllb ) 
as in the spherical case IjMeliani et al. 2 0041. It increases 
even more because of the stronger expansion of the rela- 
tivistic jet in the super- Alfvenic region. In fact, in the rel- 
ativistic solution displayed in Fig. II lb . collimation starts 
at 50 Alfven radii, while, in the non relativistic solution, 
Fig. II lb . collimation occurs only at 10 Alfven radii. 

In summary let us just point out that the mass of the 
central object and the properties of the jet (acceleration, 
morphology and energy) are not simply proportional to 
each other in the context of strong gravitational fields. 

8.3. The effect of charge separation 

As we have seen, solutions obtained with this model are 
essentially thermally driven winds but collimation is ei- 
ther thermal (pressure confinement) or magnetic (toroidal 
magnetic pinching). However, conversely to the non rela- 
tivistic case an extra decollimating force exists which is 
the electric force, despite that we neglect the light cylin- 
der effects. 

The strength of the electric field results from the in- 
duction term VpB^/c and it is higher for higher magnetic 
fields. As a matter of fact it gets more important when 
magnetic effects are important and when the light cylin- 
der is closer to the streamlines. In relativistic flows where 
the poloidal velocity is of the order of the light velocity, 
the contribution in the transverse direction of this force 
increases in the super- Alfvenic domain. 

Therefore, the contribution of this force, in relativis- 
tic jets from EMRs, is of the order of the pinching force 
and the centrifugal force. Conversely, in the non rela- 
tivistic limit, the poloidal velocity remains largely sub- 
relativistic, V p <C c, and the electric field remains weak 
(E ~ VpB^/c 2 /zu — ► 0; see for comparison Figs. EJj and 
Hk). The electric field does not affect the collimation of 
the non relativistic jet and the charge separation can be 
neglected. 

Conversely, the effects of the electric force are negligi- 
ble in pressure confined jets from IMRs as the magnetic 
effects themselves are very weak or completely negligible, 
B%jw < P/w => P E<t: Pjw. 

Parallely, the effects of the electric force on the ac- 
celeration of the jet are very weak even for EMRs. The 
correction introduced on the asymptotic speed is negligi- 
ble. This is because the electric force is perpendicular to 
the streamlines and, so, it affects mainly the morphology 
of the jet. 

9. Conclusion 

In this paper, we have investigated the problem of the 
formation and collimation of relativistic jets. We have 
explored these problems by means of a semi-analytical 
model, which is the first meridional self-similar model of 
relativistic jets, including general relativistic effects. We 
constructed it on the basis of the classical model devel- 
oped in ST94 to study jets from young stars. We have 



made an extension of this model for a black hole character- 
ized by weak angular momentum, a = JG/(Ml/c) <C 1. 
Therefore, we treated the problem of GRMHD outflows 
in a Schwarzschild metric. Moreover, we concentrated our 
efforts on modelling the jet close to its polar axis. In the 
construction of this model, we were limited to describe jets 
possessing a weak rotation velocity compared to the speed 
of light and we neglected the effects of the light cylinder. 
Thus, we restricted our study to thermally accelerated jets 
with a weak contribution of the Poynting flux. However, 
in the collimation of the jet, both electric and magnetic 
terms are comparable to the inertial and pressure gradi- 
ent ones. We have also studied the collimation effects by 
magnetic and thermal forces and the decollimation effect 
of centrifugal and electric forces. Our model is restricted 
to outflow solutions only, and we do not address to the 
problem of the origin of the hot coronal plasma. 

We found that the influence of the electric force and 
the charge separation in the jet depends on the collimation 
regime. In the case of EMRs where jets are magnetically 
collimated, the electric decollimating force plays an impor- 
tant role. This force is of the order of the magnitude of the 
centrifugal force and of the magnetic pinching. Therefore, 
relativistic jets from EMRs are less collimated than their 
non relativistic counterparts. Conversely, jets from IMRs, 
where collimation is mainly of thermal origin, are not very 
sensitive to the electric decollimation. In this type of jets, 
the contribution of p e E is balanced by the increase of the 
external pressure. 

We have also used the model to compare classical jets 
to relativistic jets. We undertook this comparative study 
by changing the free parameter \x. In fact, this new pa- 
rameter gives in the relativistic model the space curvature 
which is induced by gravity near the central object. We 
used typical values of /i from /i ~ 10 -5 for Jets from Young 
Stellar Objects to \i ~ 0.1 for jets from compact sources. 
We found that the difference between these two types of 
jets is only a scaling effect on fi for [i < 10~ 2 . In this case, 
the spatial dimensions and energies are linear functions of 
\i. For fi > 10~ 2 , the relativistic effects increase together 
with the thermal acceleration of the jet. Simultaneously, 
strongly relativistic effects tend to decollimate the jets be- 
cause of the decrease of the electric current density. 

To conclude we have seen that relativistic effects and 
particularly relativistic gravity tend to enhance the ther- 
mal acceleration (as in Meliani et al. 2004) and reduce 
the magnetic collimation of the jet (as in Bogovalov & 
Tsinganos 2005). Still collimation can be obtained either 
by thermal or magnetic means but relativistic effects lower 
the efficiency of the magnetic rotator. This means that 
despite quantitative changes, we can use this generalized 
model to verify, for the classification of AGN jets, the con- 
jectures given in Sauty et al. (2001) using the non rela- 
tivistic model. Mainly we see that the observed types of 
jets from radio loud galaxies can be connected to the ef- 
ficiency of their magnetic rotator and to the environment 
of the host galaxy. By using the present model, a more 
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precise and quantitative analysis of the observed jets will 
be presented in a following paper. 
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Appendix A: Ordinary differential equations 
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Fig. 11. Plot of the energies normalized by the parameter 
[i. In (a) is shown the relativistic solution for an EMR, 
previously displayed in Fig. Et* with /i = 10 _1 . In (c) is 
shown the corresponding non relativistic solution with fi = 
10 -5 . In (b) we plot an intermediate solution with ^ — 
10~ 2 . Plots (b) and (c) are identical which shows that 
for such small values of the energies vary linearly with 
gravity while this is not true in (a). In (d) we plot the ratio 
between the asymptotic velocity and the escape velocity 
at the surface of the corona as a function of log (1///). 
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Appendix B: The forces on the plasma 

The momentum equation can be written as 

- i2 n /2{ U.V)U-VP + ^-ElE + [V ^ h h B)] xB- 7 Wlnft = 0, 
^ 2 A w/c 2 47r 4nh 

where the generalized velocity U is given by 
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In the following we give the expressions of the various terms. 
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B.4. Magnetic force 

B.4.1. Magnetic hoop stress 
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